Creating universes with thick walls 
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We study the dynamics of a spherically symmetric false vacuum bubble embedded in a true 
vacuum region separated by a "thick wall" , which is generated by a scalar field in a quartic potential. 
We study the "Farhi-Guth-Guven" (FGG) quantum tunneling process by constructing numerical 
solutions relevant to this process. The ADM mass of the spacetime is calculated, and we show that 
there is a lower bound that is a significant fraction of the scalar field mass. We argue that the zero 
mass solutions used to by some to argue against the physicality of the FGG process are artifacts of 
the thin wall approximation used in earlier work. We argue that the zero mass solutions should not 
be used to question the viability of the FGG process. 



I. INTRODUCTION 

While our universe appears to be well described by 
ACDM cosmology and slow-roll inflation, much about the 
pre-inflationary universe remains speculative. Numerous 
models rely on quantum tunneling from some previous 
state to give an inflating universe that eventually leads 
to the universe we observe today (see for example [l[). 

We consider here the Farhi-Guth-Guven" (FGG) pro- 
cess which was originally studied in the "thin wall" limit 
0,0. In this process, a bubble of false vacuum, known 
as the seed bubble, is separated by a thin domain wall 
from a region of true vacuum. Einstein's equation implies 
two distinct solutions for the motion of the bubble wall; 
the first eventually collapses while the second expands in- 
definitely The possibility of tunneling between these 
two states is considered. Although FGG consider the case 
where a seed somehow forms in Minkowski space, other 
cases were considered (for example in §) where the seed 
forms from Hawking radiation in dc Sitter space. Either 
way, the seed collapses into a black hole but hidden be- 
hind the black hole horizon is the expanding solution. 
The mass of this bubble, M, is the m parameter in the 
usual Schwarzschild metric, the "ADM" mass. 

FGG is known to dominate over Coleman-de Luccia 
type tunneling [6] and it has been argued that this pro- 
cess can produce inflating universes that do not origi- 
nate from classical singularities 0, y, 0] ■ Despite these 
features, many calculations that study tunneling in cos- 
mology (for example in the string theory "landscape" Q) 
ignore FGG, primarily because of various arguments that 
this process might not be physical 0, H| . In this paper 
we address one of the arguments against the physicality 
of the FGG process, one that involves taking the bubble 
mass M to zero@. 

The M — > limit of the thin wall formula leads to 
a prediction of the FGG process that the probability of 



transitioning from the seed bubble to the inflating bub- 
ble remains finite even as the mass of the seed bubble 
is taken to zero. This is the ultimate free lunch, since 
it implies our universe was possibly nucleated from zero 
matter Minkowski space. However, here we argue that 
this limit is an artifact of the thin wall approximation 
which assumes that the thickness of the domain wall is 
small compared to the radius of the bubble. Indeed, as 
the radius of the bubble is taken to zero (as it is in the 
M — > limit), one should expect the thin wall approxi- 
mation to breakdown. 

In this article we examine bubbles of false vacuum sep- 
arated by a "thick wall", i.e. scalar field solutions that 
interpolate between regions of true and false vacuum. We 
construct numerical solutions for the scalar field coupled 
to gravity that are relevant to the FGG process. Prob- 
ably our most important point is an extremely simple 
one: For a fixed potential the types of possible bubbles 
are limited and the M — >• cannot even be taken. So if 
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FIG. 1. A sketch of a bubble solution. 
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one has a particular scalar field potential in mind one is 
unlikely to encounter the issues raised in [6( about FGG. 

In this work we go beyond this simple point by explor- 
ing the parameter space of a general quartic scalar field 
potential (with an overall scale fixed). We find solutions 
in such potentials cannot approach the step- function type 
solutions for cj) that are assumed in the thin wall case, 
even when the potential is made as "thin wall" as possi- 
ble. Instead, the scalar field inevitably "spills over" and 
zero mass solutions are unattainable. By comparison, 
the thin wall M — > limit relies on exact Schwarzschild 
space outside the bubble while taking the bubble radius 
to zero. While it may be possible to find an exotic po- 
tential with M arbitrarily small, we show that no quartic 
potential with a fixed overall scale admits such solutions. 



II. THE THIN WALL REVISITED 
A. The setup 

Imagine embedding a spherically symmetric bubble of 
false vacuum, the seed bubble, in a region of true vacuum 
separated by a domain wall of negligible thickness with 
surface energy density a as pictured in Fig. [TJ The thin 
wall approximation assumes that the false vacuum is de 
Sitter space, the true vacuum is Schwarzschild, and the 
stress energy tensor is discontinuous at the domain wall 

The classical solutions are discussed extensively in Q 
and it suffices to repeat a few key results. The mass M 
of the bubble is the usual Schwarzschild parameter to in 
the static foliation. This mass can be rewritten as 



M = + 47rcrrVl + r 2 - A 2 r 2 - 8n 2 Ga 2 r 3 (1) 

where A is the cosmological constant, G is Newton's con- 
stant, a is the surface energy density of the wall, and 
r is the radial coordinate in the static de Sitter and 




0.5 1.0 1.5 2.0 r 

FIG. 2. A particular (f> 4 potential with true and false vacuum 
regions labeled. 



Schwarzschild foliations. Note that the naive limit r — >• 
appears to give zero mass solutions. 

One can use the junction formalism developed in [Toj . 
The basic strategy is to place a coordinate system on the 
wall and demand continuity of the metric tensor. Then 
utilizing Einstein's equation, the rescaled radial coordi- 
nate, z, obeys the following equation 

i 1 + V(z) = E 

which is identical to that of a particle moving in a one 
dimensional potential. We know that if V > E, then 
two solutions exist, but classically the particle cannot 
move across the barrier. We can, however, have quantum 
tunneling between the two solutions. 

B. The two solutions 

As previously mentioned, two possible solutions ex- 
ist for the classical motion provided M < M cr where 
M cr ~ Af 3 is the characteristic mass of the problem Q ■ 
Here f is the radius of the bubble wall. Type (a) solutions 
are bounded solutions that begin at f = 0, expand to 
some f = r max before collapsing back to zero. These so- 
lutions avoid a classical singularity, as discussed in FGG, 
because the trajectory on the Kruskal diagram crosses to 
the right of the origin and a closed "anti-trapped" sur- 
face no longer exists. This point is further elucidated in 
0. 

Second, there are bounce solutions in which r ap- 
proaches infinity in the asymptotic past, falls to some 
minimum value, and expands again to approach infin- 
ity in the asymptotic future. The Penrose theorem im- 
plies that this space-time must have emerged from an 
initial singularity, since the bubble radius grows beyond 
(A)~i. The way to avoid this classical singularity yet 
still produce an inflating universe is to consider tunnel- 
ing between the two solutions. The two solutions are of 
identical mass and thus identical energy. This is the FGG 
process. 

The tunneling probability can be calculated using a 
functional integral [2j or a canonical quantization [3[. In 
either case, the probability of tunneling between the two 
solutions remains finite as the mass of the seed bubble, 
an input parameter, is taken to zero. 

III. VISITING THE THICK WALL 

Consider a scalar field minimally coupled to gravity in 
a quartic potential, described by the following action 

S = \m 2 P J d 4 x^(R - V a 4>V b 4>g ab - 2V((f>)) (2) 

2 

where V{<j>) — Atfi 4 — ^ycf) 3 + ^<fi 2 , rrii is the inflaton mass, 
and top is the reduced Planck mass. A particular poten- 
tial is shown in Fig. [3J 
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We work in a +2 metric signature, in reduced Planck 

units where K = 1, c = 1 and mp = VSnG . One can 
then nondimensionalize the problem by rescaling the co- 
ordinates, for example by using r* = rmp. This rescales 
the potential to 



of r. At large r, demanding true vacuum implies cp — > 
since V((f> = 0) = 0. These conditions are labeled in Fig. 

m 



A. The turning point 



V{<j>) = A0 4 - 7 3 



2m 2 p 



Here cj), A and 7 are all dimensionless. In what follows, 
all coordinates and quantities are dimensionless with 
m* = mp. Keeping mi fixed allows us to explore the 
properties of the bubble solutions without allowing the 
overall scale of the potential to vanish (in that case one 
does expect solutions with M approaching zero to be pos- 
sible). Fixing mi to the value mp is convenient for the 
dynamic range of our numerical work is also a common 
choice in inflationary models. We use standard spherical, 
(t, r, 9, 4>), coordinates. Under the assumption of spheri- 
cal symmetry, the spacetime line element takes the form 



ds 2 



-a 2 (r, t)dt 2 + a 2 (r, t)dr 2 + r 2 dfl 2 



(3) 



Note that we have not forced the metric in any region 
to take the de Sitter or Schwarzschild form, although we 
do require that the spacetime is asymptotically flat at 
large r. Here r is both a coordinate and the measure of 
proper area. The stress-energy tensor of a scalar field in 
a potential is 

Tab = d a <t>d b (t> - X -g ah {d c 4>d d ^g cd + 2V(<f>)) (4) 

Defining mass in the thick- wall case is more involved, 
since we no longer have a region of exact Schwarzschild 
space or a fixed wall position where one can place an ob- 
server. Instead, we focus on the ADM mass, which is de- 
fined at spatial infinity for asymptotically flat spacetimes. 
This is the most relevant mass for tunneling calculations 
Q. This mass is defined as [HI 



M = 2?r / drr 2 
'a 




2V(4>) 



(5) 



where prime denotes differentiation with respect to r 
while dot denotes differentiation with respect to time. 

The 00 and 11 components of Einstein's equation and 
the scalar field equation arc used in the simulation, while 
the 22 equation is used as a consistency check. 

We want to find solutions for a bubble of false vacuum 
embedded in true vacuum, i.e. we want to find the radial 
profiles and time evolution of <f>, a and a. This is done 
by demanding that <p(r = r m i n ) take the value of the 
false minimum of the potential, so that V((j> m ) acts as a 
cosmological constant near the origin. We also investi- 
gated cases in which the scalar field was not initially at 
the minimum of the potential. This did not change our 
conclusions regarding M 0, which are mainly sensitive 
to the behavior of the bubble solutions as at larger values 



For the purposes of this article it suffices to examine 
the properties of the classical solutions relevant to the 
FGG process. There is no need to find the tunneling so- 
lutions and the corresponding tunneling actions to make 
our points. Furthermore, we can understand the relevant 
properties of these solutions (namely the ADM mass) 
simply by finding the solution at its turning point, which 
further simplifies our calculations. At the turning point 
a = a = <j) = but second order time derivatives are 
nonzero. In this case, equation [5] reduces to 



M = 2- I 1- 
'0 



2V(cf>) 



dr 



(6) 



Using Einstein's equation, this can be rewritten as (again 
setting TTti = mp), 



M = 



4tt 



2ra' 



1 



1 dr 



so we see that dm/dr = only for the Schwarzschild 
metric, as expected. Additionally, this formulation gives 
PvacVbubbie for the de Sitter contribution, where p vac is 
the energy density of the vacuum, which is proportional 
to A. By inspection of Eqn. [6] (which is positive definite 
since V(</>) is everywhere positive) one can see that we 
do not expect to find M = solutions, but there is no 
apparent reason why a smooth limit to zero should not 
exist. 

There is freedom to specify the spatial profile of (j> at 
the turning point, which we choose to be <j> — c/r 2 for a 
constant c. This is consistent with spherical symmetry 
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FIG. 3. Too, demonstrating the wall "thickness", at a given 
time. 
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and is sufficiently localized to maintain an asymptotically 
flat spacetime. Choosing such an ansatz simply enforces 
locality of the bubble and does not affect the generality 
of our conclusions. 



B. The two solutions 

Evolving forward in time from the turning point solu- 
tion is used to classify the solution character. The en- 
ergy density of the expanding solution expands into the 
domain as the metric functions approach de Sitter (see 
[3]) . The energy density of the collapsing solution collapses 
immediately toward the origin while the metric functions 
approach pure Schwarzschild. A plot of T 00 (t — 0), show- 
ing the "thickness" of the wall, for the expanding solution 
is given in figure [3J Plots of the field and metric func- 
tions for expanding and collapsing solutions are given in 
figures 2] and [51 respectively. 



C. Results of trying to take M — > 

Conceptually, there are two ways in which this can be 
done. From a cosmological perspective, we can fix the 
inflaton potential and attempt to take the mass to zero 
by changing the initial condition on ip. 

On the other hand, we can tune the two constants 
in the potential while keeping a fixed 4>. (As discussed 
and motivated above, we are keeping the overall scale mi 
fixed for this investigation.) We begin with a parameter 
scan over three orders of magnitude, i.e. ranging the 
values of A and 7 from 0.1 — » 10. Let the value of the 
field at the false minimum and the maximum be i'min 
and <f>max , and the potential evaluated at these points be 
Vmin and V max , respectively. Figures [6] and [Jj show how 
the mass of the collapsing turning point solution depends 
on Acj) = 4> min - 4> max and Ay = V max - V min . 



IV. DISCUSSION 

Inspecting Eqn. [5j we see that the integrand, dm/dr, 
will not be zero unless <j) is constant and V{4>) = 0, i.e. 
exactly Schwarzschild space. The crux of our argument 
is that real potentials and fields do not admit nicely sepa- 
rated solutions; the field spills over into the whole domain 
and affects the metric functions, preventing the M — > 
limit that appears to exist in the thin-wall formalism. 

We still attempted to push the mass smoothly to zero. 
However, as Figures [5] and [JJ show, we are unable to push 
the mass below about 0.1. This is with a fixed overall 
scale set by choosing mj = nip. The point is not the 
value of mi (setting (rrii/mp) 2 = 10~ 3 does not affect 
our conclusions), but that we have fixed an overall scale. 

In any spherically symmetric problem, there is the is- 
sue of what happens at r — 0. While the numerics cannot 
evolve such a point, we can make progress analytically 
by assuming we approach exact de Sitter space, in which 
the metric functions are regular at the origin. The scalar 
field potential can then be expanding about the minimum 
to second order in <fi. The problem is then analytically 
tractable and solutions give positive mass contributions. 
Thus, our calculation of M really is a lower bound. 

Lastly, there is the issue of horizons in the computa- 
tional domain. If we had exact de Sitter, there would be 
a de Sitter horizon at r — (A)~s. We explicitly chose A, 
via parameters in the potential, to be small so that the 
horizon was significantly outside of the computational 
domain. The Schwarzschild coordinate horizon is suffi- 
ciently inside r m i n for the attempted M — > calculations. 



V. CONCLUSION 

We considered classical solutions relevant to the Farhi- 
Guth-Guven tunneling process. For a generic quartic 
potential we are unable to take the mass of our turn- 
ing point solutions smoothly to zero. Other authors have 
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FIG. 4. Turning point slice of geometry and field for the 
expanding solution. Here cj> is dotted, a is dash-dotted, and 
a is dashed. 



FIG. 5. Turning point slice of the collapsing solution. Here <j> 
is dotted, a is dash-dotted, and a is dashed. 
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shown using the thin wall approximation that the FGG 
tunneling amplitude remains finite as M — > 0, and this 
strange behavior has been used to question the physical- 
ity of the FGG process. The absence of M — > solutions 
in our more realistic thick wall calculations suggest that 
the M — > behavior is an artifact of thin wall approx- 
imation and should not be used to argue that the FGG 
process is unphysical. 
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FIG. 6. Mass of the collapsing solution, evaluated at the 
turning point, as a function of A<j>. Points represent a scan 
of potential parameters A and 7. Small values of M were not 
found in the scan. 
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FIG. 7. Mass of the collapsing solution, evaluated at the 
turning point, as a function of AU. Points represent a scan 
of potential parameters A and 7. Small values of M were not 
found in the scan. 
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